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Abstract  The  relationship  between  the  Kuroshio  volume 
transport  east  of  Taiwan  (~24°N)  and  the  impinging  meso¬ 
scale  eddies  is  investigated  using  8 -year  reanalysis  of  a 
primitive  equation  ocean  model  that  assimilates  satellite 
altimetry  and  SST  data.  The  mean  and  fluctuations  of  the 
model  Kuroshio  transport  agree  well  with  the  available  ob¬ 
servations.  Analysis  of  model  dynamic  heights  and  velocity 
fields  reveals  three  dominant  eddy  modes.  The  first  mode 
describes  a  large  eddy  of  ~500  km  in  diameter,  centered  at 
-22°  N.  The  second  mode  describes  a  pair  of  the  north-south 
counter-rotating  eddies  of-400  km  in  diameter  each,  centered 
at  23°  and  20°  N,  respectively.  The  third  mode  describes  a  pair 
of  the  east-west  counter-rotating  eddies  of-300  km  in  diam¬ 
eter  each,  centered  at  21°  N.  The  associated  velocity  fields 
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indicate  eddies  extending  to  600-700  m  in  depth  with  vertical 
shears  concentrated  in  the  upper  400  m.  All  three  modes  and 
the  model  Kuroshio  transport  have  similar  dominant  time- 
scales  of  70-150  days  and  generally  are  coherent.  The  de¬ 
creased  Kuroshio  volume  transports  typically  are  associated 
with  the  impinging  cyclonic  eddies  and  the  increased  trans¬ 
ports  with  the  anticyclonic  eddies.  Selected  drifter  trajectories 
are  presented  to  illustrate  the  three  eddy  modes  and  their 
correspondence  with  the  varying  Kuroshio  transports. 

Keywords  Western  North  Pacific  •  Mesoscale  eddy  • 
Kuroshio  transport  •  EOF  •  SVD 

1  Introduction 

Mesoscale  eddies  are  abundant  in  the  ocean.  Chelton  et  al. 
(2007),  based  on  10  years  of  altimetry  sea  surface  height 
anomaly  (SSHA)  data,  indicated  that  more  than  50  %  of 
low-frequency  variability  over  much  of  the  ocean  can  be 
accounted  for  by  eddies  with  diameters  of  100-200  km  and 
amplitudes  of  5-25  cm.  These  eddies  propagate  nearly  due 
west  at  approximately  the  phase  speed  of  nondispersive 
baroclinic  Rossby  waves  (Barron  et  al.  2009).  In  the  western 
North  Pacific,  the  most  energetic  eddy  activity  is  in  the 
Kuroshio  Extension.  There  is  yet  another  zonal  band  in 
the  subtropical  gyre  between  19°  and  25°  N  where 
energetic  eddy  activities  are  attributed  to  the  baroclinic 
instabilities  in  the  Subtropical  Countercurrent  (e.g.,  Qiu 
1999;  Kobashi  and  Kawamura  2002;  Qiu  and  Chen  2010; 
Rudnick  et  al.  2011).  The  dominant  eddy  time  scale  in  the 
STCC  is  -100  days  and  their  westward  propagating  speed  is 
-8  km/day  (Qiu  and  Chen  2010).  Most  of  the  eddies  are 
concentrated  near  22°  N,  and  many  may  reach  the  Kuroshio 
off  the  eastern  coast  of  Taiwan  (Hwang  et  al.  2004). 
Roemmich  and  Gilson  (2001)  presented  hydrographic  evi¬ 
dence  for  mesoscale  eddies  in  the  STCC  based  on  repeat 
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Fig.  1  a  Mean  model  dynamic  heights  referenced  to  1,000  dbar.  The 
red  box  marks  the  ETC  shown  in  b.  b  Topographic  map  of  the  ETC 
(depth  contour  interval  within  1,000  m  is  200  m).  The  transport  inte¬ 
grated  section  is  marked  by  the  red  line.  Mooring  stations  M1-M4  of 


PCM-1  are  marked  by  the  blue  dots,  c  Longitudinal  cross-section 
(144°E,  dotted  blue  line  in  a)  of  temperature  {blue  contours  in  °C) 
and  zonal  velocity  {black  contours  in  m/s,  positive  is  eastward) 
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Fig.  1  (continued) 


high-resolution  expendable  bathythermograph  (XBT)  tran¬ 
sects  between  Guam  and  Hong  Kong.  They  found  an  eddy- 
rich  latitude  at  -20°  N.  The  typical  eddy  size  is  about  500  km 
with  a  peak-to-trough  temperature  difference  of  2.2  °C  in  the 
center  of  the  thermocline.  Lee  et  al.  (2003)  investigated  the 
time  history  and  vertical  structure  of  a  large  cyclonic  eddy 
using  XBT  and  Topex/Poseidon  altimetry  data.  They  sug¬ 
gested  that  cyclonic  eddies  could  be  also  generated  by  slowly 
moving  typhoons. 

The  impact  of  eddies  on  the  Kuroshio  volume  transport  has 
long  been  recognized  (Yang  et  al.  1999;  Zhang  et  al.  2001). 
The  Kuroshio  originates  east  of  the  Philippines,  flows  along 
the  east  coasts  of  Luzon  and  Taiwan,  and  enters  the  East  China 
Sea  (ECS)  through  the  East  Taiwan  Channel  (ETC),  a  narrow 
(<150  km)  ridge  between  the  northeastern  coast  of  Taiwan 
and  the  southern  tip  of  the  Ryukyu  Islands  (Fig.  1).  During  the 
World  Ocean  Circulation  Experiment,  the  Kuroshio  volume 
transport  was  monitored  in  the  ETC  with  a  dense  array  of  current 
meter  moorings,  the  PCM-1  (Fig.  lb),  for  about  20  months  from 
September  1994  to  May  1996  (Johns  et  al.  2001;  Zhang  et  al. 
2001).  The  observed  volume  transport  had  a  mean  of  21.5  Sv 
(1  Sv=10^  w?/s)  and  signifrcant  100-day  timescale  fluctuations 
with  amplitudes  as  large  as  10  Sv.  For  comparison,  the  Florida 


Current  has  a  larger  mean  transport  (31.7  Sv)  but  much  smaller 
eddy  timescale  fluctuations  (Leaman  et  al.  1995;  Lee  et  al.  2001). 
The  difference  between  these  two  western  boundary  current 
systems  suggests  that  the  Kuroshio  system  is  strongly  impacted 
by  the  westward  propagating  mesoscale  eddies  (Yang  et  al.  1999; 
Zhang  et  al.  2001;  Gawarkiewicz  et  al.  2011). 

Qiu  and  Chen  (2010)  described  the  “eddy  rich”  and  “eddy 
weak”  years  in  the  STCC  (19-25°  N,  135-175°  E).  Chang 
and  Oey  (2011)  found  that  the  interannual  Kuroshio  trans¬ 
ports  estimated  from  sea  level  differences  between  Keelung 
and  Ishigaki  are  highly  correlated  with  the  STCC  variability. 
The  Kuroshio  transport  lags  the  STCC  variation  by  6- 
8  months  at  -132°  E,  and  the  time  lag  increases  eastward. 
Chang  and  Oey  (2012)  introduced  the  Philippines-Taiwan 
Oscillation  (PTO),  a  climate  index  that  relates  the  interannual 
variation  of  the  STCC  eddy  activities  and  Kuroshio  trans¬ 
ports  to  the  large-scale  wind  stress  curl  over  the  subtropical 
western  North  Pacific  (also  see  Qiu  and  Chen  2010).  The 
Kuroshio  Current  also  experiences  strong  interannual  vari¬ 
ability  further  downstream  in  the  ECS  (Liu  and  Gan  2012). 
Ichikawa  et  al.  (2008)  and  Hsin  et  al.  (2011)  showed  that  the 
Kuroshio  variations  in  ECS  were  induced  by  mesoscale 
eddies  east  of  Taiwan. 
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The  relationship  between  the  Kuroshio  transport  and  me- 
soscale  eddies  east  of  Taiwan  was  examined  from  PCM-1 
observations.  Yang  et  al.  (1999)  suggested  that  impinging 
anticyclonic  eddies  resulted  in  an  increase  of  the  Kuroshio 
transport.  Zhang  et  al.  (2001),  on  the  other  hand,  found  that 
anticyclonic  eddies  caused  a  decrease  in  the  Kuroshio  trans¬ 
port.  Also,  Gawarkiewicz  et  al.  (2011)  found  a  maximum 
correlation  of  0.83  between  PCM-1  Kuroshio  transport  and 
SSHA  at  23.9°  N,  123.2°  E.  There  were  other  studies  that 
utilized  direct  velocity/transport  measurements.  Andres  et  al. 
(2008a)  illustrated  three  events  of  the  approaching  anticy¬ 
clonic  (cyclonic)  eddies  that  increased  (decreased)  the  trans¬ 
ports  in  the  ECS  through  the  Kerama  Gap.  Also,  Zhu  et  al. 
(2004)  showed  that  the  approaching  anticy clonic/cyclonic 
eddies  increased/decreased  the  volume  transport  southeast 
of  Okinawa  Island. 

Maps  of  sea  level  anomaly  such  as  those  produced  by 
Archiving,  Validation  and  Interpretation  of  Satellite  Oceano¬ 
graphic  data  (AVISO)  are  often  used  in  studies  of  the  large- 
scale  flow  variability  (e.g.  Centurioni  et  al.  2008;  Maximenko 
et  al.  2009).  The  altimetry  spatial  resolution,  however,  is  too 
coarse  to  map  the  Kuroshio  transport.  During  the  Jason  1- 
Topex/Poseidon  tandem  mission,  the  spatial  resolution  is  im¬ 
proved,  but  the  record  is  too  short  (Scharffenberg  and  Stam¬ 
mer  2010).  The  direct  Kuroshio  volume  transport  measure¬ 
ment  is  difficult,  and  the  record  length  is  too  short  compared  to 
the  typical  eddy  timescales  (Johns  et  al.  2001).  In  this  study, 
we  analyze  the  Kuroshio  transport  variability  using  8  years 
(2003-2010)  of  reanalysis  product  from  a  high-resolution 
regional  ocean  model,  the  East  Asian  Seas  Nowcast/Forecast 
System  (EASNFS)  of  the  US  Naval  Research  Laboratory 
(NRL)  (Rhodes  et  al.  2002;  Ko  et  al.  2009)  (http:// 
www7320.nrlssc.navy.mil/NLIWI_WWW/EASNFS_ 
WWW/EASNFS_intro.html).  Unlike  the  previous  ~20-month 


PCM-1  transport  measurements,  the  8-year  EASNFS 
reanalysis  time  period  is  considerably  longer  than  the  eddy 
timescales.  Our  goal  is  to  identify  patterns  of  eddy  vari¬ 
ability  east  of  Taiwan  to  better  describe  their  relationship 
with  the  Kuroshio  transport  variability.  We  demonstrate 
that  the  mean  and  fluctuations  of  the  Kuroshio  transport 
from  the  reanalysis  are  fully  consistent  with  the  statistics 
derived  from  the  PCM-1.  We  show  that  the  eddy  kinetic 
energy  (EKE)  from  the  reanalysis  agrees  with  that  calcu¬ 
lated  from  the  altimetry.  Moreover,  we  validate  the  model 
eddy  flow  fields  with  selected  drifter  trajectories  from  the 
Global  Drifter  Program  (GDP).  Unlike  previous  studies 
based  on  either  sea  level  anomalies  or  short-volume  trans¬ 
port  measurements,  our  study  provides  a  much  more  co¬ 
herent  description  of  the  Kuroshio  transport  and  the  im¬ 
pinging  mesoscale  eddies  east  of  Taiwan.  In  particular,  we 
are  able  to  resolve  the  long-standing  controversy  from  the 
PCM-1  study  about  the  role  of  the  anticyclonic/cyclonic 
eddies  on  the  Kuroshio  transport  variations. 

2  Methods 

The  EASNFS  is  an  application  of  NRL  Ocean  Nowcast/ 
Forecast  System  (Ko  et  al.  2008).  The  ocean  model  in  the 
EASNFS  is  the  Navy  Coastal  Model  (Martin  2000)  which  is 
adopted  from  the  Princeton  Ocean  Model  (POM)  with  mod¬ 
ifications  to  accommodate  data  assimilation,  hybrid  vertical 
coordinates,  and  multiple  nesting.  A  statistical  regression 
model,  the  Modular  Ocean  Data  Assimilation  System  (Carnes 
et  al.  1996;  Fox  et  al.  2002),  based  on  historical  observations, 
is  used  to  produce  the  three-dimensional  ocean  temperature 
and  salinity  analyses  from  satellite  altimetry  and  Multi- 
Channel  Sea  Surface  Temperature  (MCSST).  EASNFS 
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Fig.  2  Mean  EKE  derived  from  model  results  (left)  and  AVISO  altimetry  (right).  Contour  lines  range  from  0  to  600  with  an  interval  of  50.  The  unit  is 
em^/s^.  Values  greater  than  800  are  assigned  to  the  same  eolor  of  800 
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assimilates  the  analyses  by  continuous  modification  of  model 
temperature  and  salinity  toward  the  analyses  using  a  vertical 
weighting  function  that  reflects  the  oceanic  temporal  and 
spatial  correlation  scales  and  the  relative  confidence  between 
the  model  forecasts  and  analyses  (Chapman  et  al.  2004;  Ko 
et  al.  2008). 

Surface  forcing  of  wind  stress,  heat  flux,  solar  radiation,  and 
surface  atmospheric  pressure,  is  derived  from  the  Navy  Opera¬ 
tional  Global  Atmospheric  Prediction  System  (Rosmond  1992). 
In  the  heat  flux,  solar  radiation  is  retained  separately,  while  the 
rest  is  adjusted  by  the  differences  between  model  sea  surface 
temperature  (SST)  and  MCSST  and  between  model  SST  and 
seasonal  climatology.  The  surface  salt  flux  is  estimated  from  the 
differences  between  model  sea  surface  salinity  (SSS)  and  the 
analysis  and  between  model  SSS  and  seasonal  climatology  (Ko 
et  al.  2008).  Jerlov  (1964)  type  lA  water  is  applied  for  solar 
extinction.  Relaxation  time  of  5  and  30  days  is  used  for  adjust¬ 
ment  of  MCSST  and  seasonal  SST  climatology,  respectively,  and 


relaxation  time  of  15  and  30  days  is  used  for  adjustment  of  SSS 
and  seasonal  SSS  climatology,  respectively. 

The  full  EASNFS  modal  domain  covers  from  17.3°  S  to 
52.2°  N  and  from  99.2°  E  to  158.2°  E.  The  horizontal  resolu¬ 
tion  is  -1/12  degree  that  ranges  from  -9.8  km  at  the  equator  to 
-6.5  km  at  the  model’s  northern  boundary.  The  horizontal 
resolution  is  about  9  km  in  the  study  region.  Temporal  reso¬ 
lution  is  the  5 -day  averaged  model  field.  There  are  41  sigma-z 
levels  with  denser  levels  in  the  upper  water  column  to  better 
resolve  the  upper  ocean.  EASNFS  analysis  has  been  previ¬ 
ously  applied  for  several  studies  in  the  western  North  Pacific, 
e.g.,  the  currents  in  the  Korea/Tsushima  Strait  (Jacobs  et  al. 
2005;  Teague  et  al.  2006),  transports  in  the  Gulf  of  Papua 
(Keen  et  al.  2006),  impact  of  typhoons  on  the  upper  ocean 
temperature  (Lin  et  al.  2008),  and  thermal  structures  in  South 
China  Sea  (Chang  et  al.  20 10). We  used  the  first  8  years 
(2003-2010)  of  the  EASNFS  reanalysis  product  in  this  study. 
The  sea  surface  height  (SSH)  from  the  EASNFS  analysis  is 
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Fig.  3  Mean  (left)  and  standard  deviation  (right)  of  veloeities  (U,  V)  and  temperature  (T)  along  PCM-1  from  model  results.  PCM-1  stations  Ml  and 
M4  eorresponding  to  eaeh  subplot  are  marked  on  the  first  subplot 
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highly  correlated  (y=0.84)  with  the  altimetry  gridded  data.  In 
this  study,  the  dynamic  height  calculated  from  the  EASNFS 
temperature  and  salinity  fields  is  used  to  describe  the  SSH  in 
order  to  avoid  effects  due  to  large-scale  barotropic  fluctua¬ 
tions.  The  spatial  and  temporal  patterns  of  dynamic  heights  are 
determined  using  the  Empirical  Orthogonal  Functions  (EOF) 
analysis.  Also,  the  relationship  between  dynamic  heights  and 
model  currents  are  analyzed  with  a  covariance  analysis  that  is 
effective  in  extracting  the  coherent  structures  of  two  diverse 
fields.  The  method  was  first  introduced  into  climate  research 
by  Bretherton  et  al.  (1992),  and  the  term  “singular  value 
decomposition  (SVD)  analysis”  was  used.  The  oceanographic 
application  can  be  found,  for  example,  in  Wang  et  al.  (2003). 

3  Results 

Figure  la  shows  model  mean  dynamic  height  (reference 
level=  1,000  db)  averaged  over  8  years  (2003-2010)  in  10- 
30°  N  and  120-150°  E.  The  Kuroshio  and  the  North  Equa¬ 
torial  Current  (NEC)  are  clearly  identified;  the  maximum 


dynamic  height  is  about  230  cm,  located  east  of  Ryukyu 
Islands.  The  NEC  bifurcation  near  14°  N  is  also  quite  no¬ 
ticeable.  The  model  mean  dynamic  height  agrees  well  with 
the  climatological  estimates  of  Qiu  (1999)  and  Rio  et  al. 
(2011),  with  significant  details  such  as  the  strength  of  the 
STCC  feature  generally  lying  within  the  range  of  estimates. 
Figure  Ic  shows  a  cross-section  of  mean  temperature  and 
zonal  velocity  at  144°  E.  The  STCC  at  24°  N  is  indicated  by  a 
weak  eastward  current  in  the  upper  300  m  and  the  corre¬ 
sponding  uplifting  isotherms.  The  model  STCC  temperature 
and  velocity  distributions  also  agree  well  with  the  observa¬ 
tions  presented  in  Qiu  and  Chen  (2010).  Figure  2  shows  a 
comparison  of  the  EKE  from  the  model  and  altimetry.  Both 
are  averaged  over  8  years  (2003-2010),  with  a  120-h  filter  to 
remove  the  short-term  fluctuations.  The  two  have  similar 
patterns  except  in  the  Kuroshio  region  where  the  altimetry 
is  known  to  underestimate  the  EKE  (e.g.,  Centurioni  et  al. 
2008),  mainly  because  of  its  inadequate  spatial  resolution. 
The  spatial  resolution  of  the  merged  data  set  is  estimated  to 
be  -150  km  by  Ducet  et  al.  (2000),  which  is  about  the  same 
as  the  width  of  Kuroshio  at  the  region. 
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Fig.  4  The  Kuroshio  volume  transport  at  PCM-1  from  model  results 
(uppermost  panel,  unit  in  Sv)  and  latitude-time  plots  of  model  dynamie 
heights  for  longitudes  133“123°E.  Mean  tmneated  PCMl  transport  is 
marked  by  the  dotted  line  at  the  uppermost  panel.  The  tilted  solid/ dashed 


lines  mark  the  propagating  eddies  related  to  the  Kuroshio  lows/highs. 
Every  January  1  from  2003  to  2010  in  the  time  axis  are  marked  with 
numbers  03-10  and  also  dots  on  the  top  panel 
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The  model  output  is  interpolated  to  the  mooring  locations 
(Ml  to  M4)  of  the  PCM-1  array  (Fig.  lb;  Zhang  et  al.  2001). 
To  compare  with  observation  (Fig.  2  in  Zhang  et  al.  2001), 
Fig.  3  shows  across-stream  structures  of  model  mean  and 
standard  deviation  for  eastward  (U)  and  northward  (V)  ve¬ 
locities  and  temperature  (T).  The  mean  V  field  shows  the  core 
(“main  axis”)  of  the  Kuroshio  about  100  km  in  width  and 
400  m  in  depth,  located  in  the  western  half  of  the  channel. 
The  mean  U  field  is  dominated  by  a  weak  eastward  flow  over 


the  entire  section,  with  larger  values  (~  0.3  m/s)  in  the 
western  half  of  the  channel.  The  base  of  the  main  thermo- 
cline  (16  °C)  rises  from  about  300  m  in  depth  at  the  eastern 
side  to  200  m  at  the  western  side.  The  corresponding  vertical 
shears  are  also  stronger  on  the  western  side.  Standard  devi¬ 
ations  of  U  and  V  are  comparable  along  the  entire  section. 
The  standard  deviation  of  U  increases  eastward  and  inten¬ 
sifies  near  the  surface,  and  the  standard  deviation  of  V  has  a 
core  (-0.2  m/s)  in  the  western  half  of  the  channel.  The 


a 


Fig.  5  a  EOF  amplitudes  for  dynamic  heights:  first  mode  (left),  second 
mode  (middle),  and  third  mode  (right).  EOF  amplitude  is  normalized 
and  the  negative  region  is  shaded.  Percentages  of  each  mode  are  also 
shown,  b  EOF  principal  components  (EOF_PCl-3)  for  dynamic  heights 


(top  three  plots)  and  Kuroshio  transport  (black,  bottom  plot,  also  shown 
in  Fig.  4).  Three  vertical  dashed  lines  marked  at  days  951,  2,331,  and 
2,851  correspond  to  modes  1,3,  and  2  in  Fig.  10 
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standard  deviation  of  T  is  large  in  the  surface  layer,  increas¬ 
ing  eastward.  The  model  velocity  and  temperature  structures 
agree  well  with  the  PCM-1  observations  (Zhang  et  al.  2001). 

The  model’s  Kuroshio  volume  transport  through  the  ETC 
is  integrated  along  the  red  line  marked  in  Fig.  lb.  It  ranges 
from  10  to  39  Sv  ,with  dominant  timescales  of  70-150  days 
(Fig.  4).  The  8-year  mean  volume  transport  is  24.4  Sv,  with  a 
standard  deviation  of  4.7  Sv.  For  the  truncated  PCMl  section 
(M1-M4),  the  model  transport  is  2 1.8 ±5. 3  Sv,  which  agrees 
well  with  2 1.5  ±4.1  Sv  obtained  from  the  20-month  PCM-1 
observations  (Johns  et  al.  2001).  We  noted  that  the  PCM-1 
period  is  relatively  short  and  is  not  included  in  the  reanalysis. 
The  good  agreement  suggests  that  the  Kuroshio  transport  is 
quite  stationary  in  the  interannual  time  scale.  The  total  mean 
transport  of  24.4  Sv  also  agrees  with  the  Kuroshio  transport 
in  the  East  China  Sea  at  28°  N,  24.0±0.9  Sv,  based  on  13- 
month-long  (November  2003  to  November  2004)  combined 
velocity  and  sea  level  observations  (Andres  et  al.  2008b). 
Figure  4  also  displays  the  latitude-time  plot  of  the  8-year 
model  dynamic  heights.  To  include  the  shallower  region,  the 
reference  level  for  dynamic  heights  is  reduced  to  500  m 
(DH_500).  The  boxes  from  the  bottom  to  the  top  display 
the  dynamic  height  contours  in  17-25°  N  and  123-133°  E. 
Westward-propagating  eddies  can  be  clearly  identified  from 
the  time  displacement  in  longitude.  Zhang  et  al.  (2001) 
suggested  that  the  low  Kuroshio  volume  transports  are  asso¬ 
ciated  with  the  approaching  anticyclonic  eddies.  In  Fig.  4, 
several  transport  highs  and  lows  are  traced  back  to  the 
approaching  eddies.  Our  results  indicate  that  the  low  (high) 
Kuroshio  volume  transports  generally  are  associated  with  the 
westward-propagating  cyclonic  (anticyclonic)  eddies.  Only  a 
few  cases  agree  with  Zhang  et  al.  (2001),  for  example,  the 
last  case  illustrated  in  Fig.  4  in  which  an  approaching  anti¬ 
cyclonic  eddy  caused  a  decreased  Kuroshio  transport. 


Fig.  6  Normalized  energy  speetra  of  EOF_PCl-3  and  model  Kuroshio 
transport 


The  dynamic  height  has  a  strong  annual  signal  due  to  steric 
variations  (Wang  andKoblinsky  1996;  Stammer  1997).  Since 
our  study  is  primarily  focused  on  mesoscale  eddies,  the  annual 
harmonics  which  accounted  for  40  %  of  the  total  variance  are 
first  removed.  EOF  analysis  is  then  applied  to  the  residual 
dynamic  heights  (rDH_500)  in  the  model  sub-domain  bound¬ 
ed  byl9-24°  N  and  122.7-127.2°  E.  Figure  5  shows  spatial 
(amplitudes)  and  temporal  (principal  components,  PCs)  struc¬ 
tures  of  the  first  three  EOF  modes  which  explain  about  30, 20, 
and  14  %  of  the  total  variance,  respectively.  The  first  mode  is  a 
large  single  eddy  of  about  500  km  in  diameter,  centered  at  22° 
N.  The  second  mode  indicates  two  north-south  counter¬ 
rotating  eddies  of  about  400  km  in  diameter  each,  centered 
at  20°  and  23°  N,  respectively.  The  third  mode  shows  two 
east-west  counter-rotating  eddies  of  about  300  km  in  diameter 
each,  centered  at  123.5°  and  126.5°  E  of  21°  N,  respectively. 
We  noted  that  since  the  eddies  are  propagating  from  the  open 
ocean,  the  EOF  modal  structures  will  be  affected  by  the 
analysis  domain.  The  box  described  earlier  is  chosen  to  focus 
on  the  region  east  of  Taiwan  and  upstream  of  ETC,  as  only 
when  eddies  are  close  to  ETC  could  high  coherence  with  the 
Kuroshio  transport  be  expected  (Gawarkiewicz  et  al.  2011; 
Chang  and  Oey  2011).  The  large-scale  eddy  structures  in 
STCC  are  examined  later  (Fig.  8). 

Figure  6  shows  the  spectra  of  the  first  three  principal 
components  (PCI,  PC2,  and  PC3;  PCl-3)  and  the  Kuroshio 
transport  in  the  ETC.  All  four  spectra  show  dominant  energy 
peaks  at  around  70-150  days.  Figure  7  displays  the  coherency 
and  phase  lags  of  the  Kuroshio  transport  with  PC  1-3.  For  the 


Frequency  (1 /day) 

Fig.  7  Coherency  and  phase  between  EOF  PC 1  -3  and  model  Kuroshio 
transport 
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Fig.  8  Correlation  coefficients  between  dynamic  heights  and  EOF_PCl-3  from  -120-day  to  +30-day  lags  with  an  interval  of  30  days:  first  mode 
(left),  second  mode  (middle),  and  third  mode  (right) 
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Fig.  9  Modal  structures  of 
dynamic  height  with  three- 
dimensional  velocity  field  from 
SVD  analysis.  Left,  dynamic 
height  {color  shading  with 
contour  lines)  and  surface 
velocity  (vectors)  of  the  first 
three  modes.  Right,  cross-section 
of  normal  velocity  along  the  lines 
marked  in  the  left  panel.  Triangle 
and  multiplication  symbol  denote 
the  eddy  centers 
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first  mode  (Kuroshio  transport  vs.  PCI),  the  coherency  is 
high  for  periods  >75  days,  with  the  transport  lagging  PCI 
by  ~20  days.  There  is  also  a  high  coherency  at  40  days 
with  the  transport  leading  PCI  by  ~10  days.  For  the  second 
mode,  the  coherency  generally  is  low  except  at  a  very  long 
period  >150  days,  with  the  transport  lagging  by  20-30  days. 
For  the  third  mode,  there  is  a  high  coherency  at  75  days 
(transport  leading  by  ~20  days)  and  40  days  (~  zero  lag). 


respectively.  The  coherency  is  also  high  at  a  very  long 
period  >150  days.  We  noted  that  although  the  EOF  modes 
are  orthogonal  in  time  domain,  they  could  be  coherent  in 
certain  frequency  bands.  For  example,  PC2  and  PCS  are 
coherent  at  70-day  period  (not  shown).  We  choose  the  time 
domain  EOF  analysis  in  line  with  the  previous  studies.  The 
time  domain  EOF  analysis  also  allows  for  easy  reconstmc- 
tion  of  the  ocean  state  (for  example.  Fig.  10). 


^  springer 


Ocean  Dynamics  (2013)  63:1027-1040 


1037 


In  the  EOF  analysis,  it  is  assumed  that  eddies  are  quasi¬ 
stationary  within  the  analysis  domain.  To  examine  how  the 
three  EOF  modes  are  related  to  dynamic  heights  in  the  open 
ocean,  the  principal  components  are  lag-correlated,  with 
residual  dynamic  heights  (rDH_500)  in  the  STCC.  Figure  8 
shows  lag  correlations  every  30  days  from  -120  days  to  + 
30  days  for  the  three  EOF  modes.  With  zero  time  lag,  the 
correlation  map  shows  spatial  patterns  that  are  basically  the 
same  as  the  EOF  modal  pattern.  This  indicates  that  the  EOF 
analysis  box  indeed  has  captured  the  structures  of  the  im¬ 
pinging  eddies.  Also,  at  zero  time  lag,  PCI  and  2  are  posi¬ 
tively  correlated  with  rDH_500  around  southeast  of  ETC 
where  Gawarkiewicz  et  al.  (2011)  found  the  maximum  cor¬ 
relation  between  the  Kuroshio  transport  and  SSHA.  With 
time  lags,  all  three  EOF  modes  can  be  tracked  over  a  long 
distance/time  in  the  STCC  between  20  and  24°  N.  The  phase 
speeds  are  ~8  km/day  (8.16,  8.03  and  7.4  km/day,  respec¬ 
tively,  for  PC  1-3).  Also,  only  the  first  mode  seems  to  have 
significant  impact  further  downstream  in  the  East  China  Sea 
(at  +30-day  lag). 

The  model’s  three-dimensional  velocity  field  is  analyzed 
for  the  modal  structure  in  the  vertical.  Instead  of  treating 
velocities  as  an  independent  field,  SVD  analysis  is  used  to 
extract  the  flow  patterns  that  are  correlated  with  the  dynamic 
height  (Bretherton  et  al.  1992).  For  consistency,  the  annual 
cycle  is  also  removed  from  the  velocity  field.  The  SVD  is 
applied  between  the  dynamic  heights  and  horizontal  veloci¬ 
ties  from  the  surface  to  1,000  m  in  depth.  Figure  9  shows  the 
surface  velocity  and  dynamic  height  patterns  for  the  first 
three  modes  and  the  corresponding  cross-sections  of  the 
normal  velocity  field.  The  time  series  of  the  first  SVD  mode 
is  similar  to  EOF_PC2  (correlation  coefficient,  y=0.78),  the 


second  SVD  mode  to  EOF_PCl  (y=0.77),  and  the  third 
SVD  mode  to  EOF_PC3  (y =0.89).  The  velocities  are  highly 
coherent  in  the  vertical.  The  flow  field  extends  to  -700  m  in 
depth  for  modes  1  and  2  and  to  -600  m  for  mode  3.  Note  that 
most  of  the  vertical  shears  are  confined  to  the  upper  400  m. 
The  deep  penetration  of  the  eddy  motion  is  consistent  with 
the  findings  of  Roemmich  and  Gilson  (2001)  who  found  that 
the  temperature  anomalies  associated  with  the  “composite” 
warm/cold  eddies  extend  to  600-800  m,  with  the  largest 
temperature  variations  (and  hence,  the  geostrophic  velocity 
shears)  in  the  upper  400  m. 

A  realization  of  the  ocean  state  is  generally  a  combination 
of  the  three  leading  EOF  modes  (64  %  of  total  variance). 
There  are  instances  though  when  a  single  mode  dominates. 
Figure  9  illustrates  examples  when  a  single  EOF  mode  is 
revealed  from  the  concurrent  Surface  Velocity  Program  drift¬ 
er  trajectories;  the  latter  were  obtained  from  the  enhanced 
GDP  database  maintained  at  the  Scripps  Institution  of 
Oceanography  (Niiler  2001).  For  each  case,  a  snapshot  of 
the  model  dynamic  height  (DH_500)  is  presented  (with  the 
area  average  removed  for  clarity),  and  all  available  drifter 
tracks  within  ±30  days  of  the  snapshot  are  superimposed. 
The  time  corresponding  to  each  snapshot  is  marked  in  the 
EOF  time  series  of  Fig.  5b.  The  first  case  is  at  day  951 
(2005/8/9)  when  the  first  EOF  mode  dominated.  The  situa¬ 
tion  is  similar  to  SVD  mode-2  (Fig.  9)  in  that  a  large  cyclonic 
eddy  (negative  PCI)  is  present  east  of  Taiwan.  Two  drifters 
(“red”  and  “blue”)  closely  traced  the  cyclonic  eddy.  The  third 
drifter  (“black”),  however,  was  carried  in  a  large  anticyclonic 
circulation  south  of  the  cyclonic  eddy.  The  second  case  is  at 
day  2851  (2010/10/22)  when  the  second  EOF  mode  domi¬ 
nated  (negative  PC2).  The  drifters  initially  were  all  trapped 


Fig.  10  Drifter  trajectories  from  the  GDP  database  and  corresponding 
modeled  dynamic  height  anomalies  at  days  951  (2005/8/9,  left),  2,851 
(2010/10/22,  middle),  and  2,331  (2009/5/20,  right)  for  demonstrating 
EOF  modes  1,  2,  and  3,  respectively.  The  domain-averaged  dynamic 
height  is  subtracted  of  each  plot.  Solid  (dashed)  lines  denote  positive 


(negative)  values  with  a  contour  interval  of  10  cm.  Trajectories  within  ± 
30  days  are  shown:  -20  days  (inverted  filled  triangle),  -10  days  (filled 
diamond),  0  day  (filled  circle),  ±10  days  (filled  square),  and  ±20  days 
(filled  triangle) 
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in  a  small  cyclonic  eddy,  but  two  of  them  later  escaped  and 
followed  a  large  anticyclonic  eddy  to  the  south.  This  period 
also  corresponds  to  a  (very)  low  Kuroshio  transport  (Fig.  5b) 
as  a  large  portion  of  the  Kuroshio  flow  is  diverted  eastward  at 
the  southern  tip  of  Taiwan.  The  third  case  is  at  day  2,331 
(2009/5/20)  when  the  third  EOF  mode  was  relatively  more 
important  (negative  PC3).  At  -10  day,  the  two  drifters  were 
closer  by.  The  outer  (“blue”)  drifter  made  a  loop  around  a 
cyclonic  eddy  but  moved  straight  to  the  south  between  +10 
and  +20  days,  indicating  that  the  cyclonic  eddy  had  evolved 
into  a  north-south  elongated  shape.  Correspondingly,  the 
inner  (“red”)  drifter  made  a  spiral  trajectory  along  an  anticy¬ 
clonic  eddy.  The  north-south  orientation  of  the  eddy  pair 
(east-west  counter-rotating  eddy  pair)  at  >  +10  days  is  con¬ 
sistent  with  the  EOF  mode  3. 


4  Discussion 

Using  the  8-year  (2003-2010)  EASNFS  analysis,  the  relation¬ 
ship  between  the  Kuroshio  volume  transport  and  mesoscale 
eddies  east  of  Taiwan  is  examined.  The  EOF  analysis  of 
model  dynamic  heights  identifies  three  dominant  modes,  a 
single  large  eddy  (PCI),  a  north-south-oriented  eddy  pair 
(PC2),  and  an  east-west-oriented  eddy  pair  (PC3).  The  corre¬ 
sponding  velocity  fields  indicate  that  the  eddies  penetrate  to  a 
depth  of  600-700  m,  with  strong  velocity  shears  in  the  upper 
400  m.  The  EOF  modes  and  the  Kuroshio  volume  transport 
have  similar  timescales  of 70-1 50  days,  and  they  generally  are 
coherent  with  a  lag  of  about  20  days. 

Yang  et  al.  (1999)  suggested  that  the  SSHA  east  of  Taiwan 
were  highly  correlated  with  the  Kuroshio  transports  estimated 
from  sea  level  differences  between  Keelung  and  Ishigaki.  Sim¬ 
ilarly,  Gawarkiewicz  et  al.  (2011)  found  maximum  correlation 
between  the  Kuroshio  transports  and  SSHA  near  Ishigaki.  Yang 
et  al.  (1999)  showed  that  an  approaching  anticyclonic/cyclonic 
eddy  increased/decreased  the  Kuroshio  transport.  Their  conclu¬ 
sion  is  supported  by  Andres  et  al.  (2008a)  who  examined  direct 
transport  measurements  during  several  eddy  events.  Zhang 
et  al.  (2001),  on  the  other  hand,  suggested  that  the  impinging 
anticyclonic  eddies  caused  a  decrease  of  the  ETC  transport. 
Ichikawa  et  al.  (2008)  obtained  a  similar  result  by  correlating 
the  surface  Kuroshio  velocity  determined  from  the  high- 
frequency  radar  measurements  with  the  altimetry  data.  Our 
results  reveal  a  more  complicated  situation  of  three  different 
eddy  modes  affecting  the  Kuroshio  transport.  Yang  et  al.  (1999) 
and  Gawarkiewicz  et  al.  (2011)  can  be  related  to  PCI  and  PC2, 
and  Zhang  et  al.  (2001)  can  be  related  to  PC3.  For  PCI  and 
PC2,  part  of  the  Kuroshio  transport  is  diverted  by  the  cyclonic 
eddies  to  the  east  near  the  southern  tip  of  Taiwan  (~21°  N) 
(Fig.  10).  For  PC3,  part  of  the  Kuroshio  transport  is  diverted  by 
the  anticyclonic  eddies  off  the  ETC  along  the  east  side  of  the 
Ryukyu  Islands  (-24°  N).  As  the  first  two  modes  are  more 


frequent,  cyclonic  eddies  play  the  dominant  role  in  decreasing 
the  Kuroshio  volume  transport  across  the  ETC  (Fig.  3). 

Ichikawa  (2001)  applied  CEOF  analysis  to  study  the  impact 
of  eddies  further  north  (22-34°  N)  using  15-month  altimetry 
data.  His  first  mode  represents  the  annual  signal  (-40  %  of 
variance),  and  the  second  mode  which  accounts  for  -24  %  of 
variance  has  a  timescale  of  -6  months  and  plays  the  most 
important  role  in  affecting  the  Kuroshio  transport  in  the  Tokara 
Strait.  His  second  mode  is  related  to  the  mesoscale  eddies  at 
23°  N,  corresponding  to  our  PCI  eddy.  The  lag  correlation 
analysis  indeed  shows  that  only  the  PCI  eddy  influences  the 
Kuroshio  transport  in  the  East  China  Sea  (Fig.  8). 

The  atmospheric  reanalysis  products  such  as  the  NCEP/- 
NCAR  Reanalysis  (www.esrl.noaa.gov/psd/data/gridded/ 
data.ncep.reanalysis.html)  are  routinely  treated  as  the 
datasets.  The  reanalysis  is  advantageous  as  it  employs  a 
state-of-the-art  analysis/forecast  system  to  systematically 
assimilate  the  observations  which  are  otherwise  discrete 
in  space  and  time.  An  oceanographic  reanalysis  product, 
on  the  other  hand,  is  rarely  available  because  globaFre- 
gional  ocean  operational  systems  are  still  few.  We  demon¬ 
strate  in  this  study  the  usefulness  of  EASNFS  analysis  in 
obtaining  credible  long-term  Kuroshio  volume  transport 
estimates.  Also,  Qiu  and  Chen  (2013)  and  Chang  and 
Oey  (2012)  recently  have  emphasized  the  importance  of 
large-scale  wind  stress  curl  in  the  STCC.  Presumably,  the 
impact  of  wind  forcing  on  the  Kuroshio  volume  transport 
could  be  also  diagnosed  from  the  EASNFS  reanalysis 
which  includes  realistic  atmospheric  forcing. 
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